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Abstract 

We investigate a new model of binary fragmentation with inhibition, driven 
by the white noise. In a broad range of fragmentation probabilities, the power- 
law spatiotemporal correlations are found to arise due to self-organized crit- 
icality (SOC). We find in the SOC phase a non-trivial power spectrum of 
the temporal sequence of the fragmentation events. The 1// behaviour is 
recovered in the irreversible, near-equilibrium part of this phase. 
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The concept of self-organized criticality (SOC) |L| has attracted much attention 
as it has been demonstrated that it leads naturally to power-law distributions in 
space and time. The ubiquity of such correlations in nature and, in particular, 1/ / 
noise in time , are expected to follow from the widespread occurence of SOC in 
most dissipative many-body systems. In the fragmentation processes, the idea of 
universality comes from the observation that many experimental data on cluster- 
size (mass) distribution can be fitted by a single law n& ~ k~ T , independently of 
interactions involved or orders of scale. Values of the exponent r are known suffi- 
ciently accurately to exclude the possibility of grouping in few universality classes 
by the value of r . There is here a striking resemblance with the fragmentation - 
inactivation binary (FIB) model ||, || , which splits the set of all sequentially frag- 
menting systems into three classes. Two of them , at the transition line and in 
the shattering phase, exhibit power-law size-distribution, but the possible values of 
the exponent r are smaller than 2 in one class (on the transition line) and are 
larger than 2 in the other one (the shattering phase) |J . Despite this qualitative 
resemblance these two classes are of fundamentally different nature as manifested by 
the properties of the cluster-size fluctuations [f|, |J . It has been shown that the FIB 
model is able to reproduce quantitatively the multifragmentation results in nuclear 
heavy- ion collisions at intermediate energies . 

In this Letter, we introduce the FIB automaton which is an alternative version 
of the FIB model, and show that it evolves into a SOC state with no characteristic 
time or length scales in a broad range fragmentation probabilities. We will show 
also that the occurence of the anomalous power spectra 1/ f 13 (1 < f3 < 2) in 
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sequential fragmentation is associated with the SOC phase of the FIB automaton. 
In particular, 1/f spectrum is obtained in the limit of weak fragmentation activity. 

The FIB automaton is defined on a one-dimensional lattice of linear size N (i = 
1, . . . , N). The zth coordinate is an integer denoting a mass of the object. To each 
lattice site we assign an integer n(i) which is to represent a multiplicity of fragments 
of mass i . The dynamical rules of the FIB automaton include the perturbation 
mechanism and the relaxation rules. Let us beginn by the relaxation rules. One 
deals with clusters characterized by some conservative quantity, that we shall call 
the cluster mass . The cluster of mass N is relaxing via an ordered and irreversible 
sequence of steps. The first step is either a binary fragmentation occuring with 
the probability pp or an inactivation, with probability 1 — pp . Once inactive, 
the cluster cannot be reactivated anymore. The successful fragmentation leads to 
two fragments, say (N) — > (j) + (N — j) , with the mass partition probability 
aFj^N-j , where a is given by the normalization condition for coefficients Fij . The 
following step is either the fragmentation of the cluster (j) , with probability pp , 
or its inactivation (j) — > (j)* , with probability 1 — pp . The relaxation process 
terminates when all clusters are inactive. 

We consider the fragmentation process which is driven by the 'white noise' (Fjj = 
Fij = const). It is assumed that the fragmentation probability pf along the 
fragmentation cascade has a fixed value between and 1, independent of the size of 
the fragmenting object, except for the monomers which are not allowed to break up 
[0 . Any event , either a fragmentation or an inactivation , occurs in the time interval 
between t and t + alt after the cluster has been activated with a probability pp . 
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For pf = const , the FIB process can be mapped into the directed percolation on 
the Cayley tree which is a mean-field percolation. The fragmentation is a branching 
process, since the propagator of fragmentations is independent of mass. Each node 
is occupied with a probability 1 and at each occupied point at time I dt one 
chooses between three possibilities : fragmentation, inactivation, and 'no event' 
with respective probabilities p 2 , (1 — p) 2 and 2p(l — p) . The fragmentation 
probability is : pp = p 2 /[p 2 + (1 — p) 2 } ■ This fragmentation is also analogous to the 
process of self-avoiding random walk, because the previously activated sites of this 
tree-like process 'repel' any subsequent reactivation. At each fragmentation, a given 
cluster is replaced by two descendants and the fragmentation multiplicity increases 
by one unit. The probability Pn(S) that a FIB branching process creates exactly 
Mq fragments is for large N : 



Pn[M = S] ' ' 



2p F 2S-1 



\ s i 



[4p F (l-p F )f ~S- s/2 exp(aS) if p F < 1/2 

(1) 



with a = log[4pjr(l — pp)] . For pp > 1/2 , the branching tree has an infinite size 
(Mo = 00). At last, pf = 1/2 is the critical point of the FIB branching process as 
in this case the branching tree barely survives and a fragmentation multiplicity for 
large iV is a power law with the exponent 3/2 . Moreover, the average value of the 
multiplicity , which is < M >= (1 —pp)/(l — 2pp) for pp < 1/2 becomes infinite 
if pf > 1/2 . The value pp — 1/2 is also a critical point of the avalanche process 
for the undirected sandpile process in the mean-field limit . Notice however, that 
the sequential fragmentation is also critical from the point of view of the shattering 



transition for any value of Pf , even though the branching tree size-distribution is 
modified. 

The FIB automaton is perturbed by adding one mass unit to the randomly cho- 
sen fragment (j) on the lattice {n(i), i = 1, ...,iV}. The cluster of mass N is 
unstable and, therefore, the fragmentation avalanche is started if the perturbation 
produces a cluster N . Once we have chosen the perturbation mechanism and the 
relaxation algorithm then the algorithm of the temporal evolution of the FIB au- 
tomaton goes as follows : 

(1) Specify an initial configuration {n(i), i — 1, . . . , iV} . 

(2) Whenever n(N) ^ , the cluster of size N breaks into smaller fragments accord- 
ing to the relaxation algorithm of FIB model until all branches of the FIB cascade 
become inactive. 

(3) The fragment multiplicity at all sites is updated, increasing by the amount of 
produced inactive clusters. 

(4) Choose a fragment (j) at random and increase its mass by one unit: j — > j + 1 . 
Change the fragment multiplicities at sites j : n(J) — > n(j) — 1 and j + 1 : 
n(j + 1) — > n(j + 1) + 1 . Return to step 2 . 

For p F < 1/2 , after a transient period, a stationary state of FIB automaton is 
reached in which the normalized average size of the cluster a(t) = (l/N) J2j j n(j) ex " 
hibits small fluctuations around its asymptotic value : 

a= lim 1 f T dta(t) = I ~ ^ F . (2) 
t^oo T Jo 2 — 2p F K J 

This is the 'high- viscosity' regime, where the accumulation mechanism of the FIB 
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automaton compensates exactly the relaxation mechanism and the system drives 
alone into a steady state with parameters independent of the lattice size. We shall 
see below that the steady state of the FIB automaton in this regime has all features 
of SOC state. In the upper part of Fig. 1 we show for p F — 0.1 the dependence 
of the average size of the clusters on time, where a unit time-step is one update 
of the whole lattice {n(i), i = 1, . . . , N} . For p F > 1/2 , the asymptotic average 
cluster-size a depends explicitely on the lattice size and equals ~ (lniV) -1 for 
p F = 1/2 and ~ l/[(2 - 2p F )({2p F )]N 1 - 2pF for p F > 1/2 . Thus, a becomes equal 
zero in the thermodynamic limit N — > oo . This is the 'low-viscosity' regime of 
the FIB relaxation mechanism, where the cluster-mass accumulation is completely 
destroyed. We have found , that for any p F , the cluster-size distribution in the 
steady state of the FIB automaton exhibits a power-law behaviour. However, in 
contrast to the original FIB model |J , it does not manifest any finite-size effects 
even in small lattices (see the lower part of Fig. 1). The average slope of the cluster- 
size distribution in the steady state exhibits small fluctuations about its asymptotic 
value f = 2p F . 

Before we discuss the statistical properties of FIB automata , let us first define 
the instanteneous dissipation rate for the avalanche e : 

/(e) (*) =£xi 6) (*) (3) 
k 

/( e ) (t) is an indicator function of unstable clusters at a time t in the fragmentation 
avalanche. xjj. (t) i n © is the characteristic function of the cluster k and equals 
either 1 , for t e [£fc,ifc] where tf~ is the instant of formation of the cluster k and 
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t' k the time when it disappears, or otherwise. The summation goes over all 
clusters in the avalanche e . At a given time t , the value of the characteristic 
function f(t) is then just the number of active clusters. 

The lifetime distribution of fragmentation avalanches are shown in Fig. 2 for few 
selected values of pp ■ One may notice a power-law behaviour P[T — t] ~ t 1 ^ 01 for 
p F < 1/2 , i.e. in the 'high- viscosity' regime. The values of the exponent a for 
different values of pp are summarized in Table 1. 

The avalanche-size is defined as the total dissipation of the fragmentation avalanche : 

s= /(*)<** = £ (t'k-tk) ~ M T , (4) 

This is just the sum of the lifetimes of all the clusters which have appeared in the 
sequence of breaks, where M and T are the average multiplicity and the total 
fragmentation time for a given event respectively. For pp < 1/2 , the value of 
M is a function of pp but does not depend on the size of the fragmenting object. 
Hence, we expect that s ~ t 71 with 7x = 1 . This is verified by the calculated 
size-distribution, which exhibits also a power-law behaviour P[S — s] ~ s 1_r with 
r — a . For pp > 1/2 , i.e. in the 'low- viscosity' regime, the branching tree of the 
fragmentation process survives until the low mass cutoff for monomers. In this case, 
the distribution of both avalanche-lifetimes and avalanche-sizes loose their scale- 
invariant features and are dominated by the cutoff scale. Consequently, only for 
relatively small time and size scales, the distributions P[T = t] and P[S — s] are 
consistent with a power-law behaviour. In this range of scales, a — r — 2 , 
independent of pp . 
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The sequential fragmentation in FIB automaton has properties of the SOC state 
for all pf in the 'high- viscosity' regime < pp < 1/2 . For each pp , the stationary 
state of the FIB automaton has different features as manifested e.g. by different 
asymptotic values of the average cluster-size a or the average slope of the cluster- 
size distribution f . For each value of the fragmentation probability pp (< 1/2), 
FIB model provides a different mean-field realization of a SOC state with its par- 
ticular exponents of the power-law spatiotemporal distributions. In this sense, it is 
appropriate to speak about a SOC phase of FIB automaton for < pp < 1/2 . The 
extreme points: pf = and 1/2 at the two ends of the SOC phase have a par- 
ticular significance as they correspond to the stationary limit of the fragmentation 
process and to the critical point of the branching process respectively. Below, we 
shall show, that they play also a special role in characterizing the anomalous power 
spectra S(f) ~ 1/f 13 with 1 < /3 < 2 . 

One of the main aims behind SOC studies was the search for a generic explanation 
of 1/f flicker noise|| whose widespread occurence is one of the great unresolved 
problems in physics. The concept of SOC state is an attempt to find a general 
principle behind ubiquity of 1/f noise in apparently unrelated physical systems. As 
pointed out before 0, [10], |TTJ , the spatiotemporal scaling in the SOC state does not 



necessarily manifest itself in the nontrivial features of the noise spectrum. Indeed, 
the power spectrum of the cellular automaton model of Bak, Tang and Wiesenfeld|l|] 
is 1/f 2 as for the random walk. 

Power spectrum of linearly superimposed fragmentation avalanches can be cal- 
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culated directly. For that let us define the total dissipation rate : 



(5) 



e 



where T( e ) is a random instant of time when a fragmentation in the event e has 
started and the summation goes over all events generated by the FIB automaton. 
The power spectrum of j(t) : 



is shown in Fig. 3 for few values of p F . The function S(f) indeed shows a power law 
behaviour ~ l//' 3 even in relatively small systems (N = 2048). The values of the 
power law exponent (3 for different values of p F are shown in Table 1. For p F > 1/2 , 
the power spectrum is those of a random walk, i.e. 1/ f 2 . For p F < 1/2 , i.e. in the 
SOC phase of the FIB automaton, the power spectrum is anomalous (1 < (5 < 2 ) 
in the frequency range /o < / < A ■ For / < f , one finds a 'white noise' power 
spectrum and the low-frequency limit fo is independent of N . For f > fi , the 
power spectrum becomes 1/ f 2 and the high-frequency limit f\ grows linearly with 
N . High frequencies correspond to short fragmentation times, where one is sensible 
to the fragmentation of the first cluster of size iV . Since this process is unique, 
one expects 1/ f 2 spectrum as soon as / 3> 1/t n , where t n ~ Pf/{N — 1) is a 
typical time of disappearing of the first cluster in the fragmentation avalanche. The 
exponent j3 becomes 1 in the limit p F — > , i.e. when approaching the stationary 
limit of the FIB automaton. The exponent (3 of the power spectrum S(f) is related 
to the exponents a , r of the lifetime- and size- distributions by the relation : 




2 



2 



> 




= (4-r) 



2-r' 



(7) 
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This relation has been derived by Christensen et al. (see eqs. (36), (39) and (40) 
of ref . [|1(| ) for sandpile cellular automata. Numerical verification of this relation is 
shown in Table 1. The validity of (0) for pp > 1/2 is less certain as the power-law 
behavior in distributions P[T = t] and P[S = s] is here restricted to relatively 
short times only. 

The proposed mechanism for anomalous activity spectrum due to critical frag- 
mentation is a 'bulk' phenomenon associated with the 'mass' fragmentation. It is an 
irreversible phenomenon, associated with the dynamical SOC phase of the fragmen- 
tation process. The coherence in the activity spectrum over all time scales (/3 — > 1) 
is found in the limit of dying out fragmentation activity (pp —>■ 0), i.e. close to 
the stable equilibrium of the system without external driving. The ubiquity of the 
1/ / phenomena follows once it has been established that the near-equilibrium dy- 
namics of most dissipative many-body systems is determined by the SOC state which 
is an attractor to the complicated many-body dynamics. In contrast to the expla- 
nation of the flicker noise, based on the deterministic spring - block models[|T2| or 



the lattice-gas model |I3| the present model has also an inherent local conservation 



law, the cluster 'mass' in each fragmentation vertex. 



10 



References 

[1] P. Bak, C. Tang and K. Wiesenfeld, Phys. Rev. Lett. 59, 381 (1987); Phys. 
Rev. A38, 364 (1988); 

C. Tang and P. Bak, Phys. Rev. Lett. 60, 2347 (1988); 

P. Bak, K. Chen and C. Tang, Phys. Lett. A147, 297 (1990). 

[2] P. Dutta and P.M. Horn, Rev. of Modern Physics 53, 497 (1981); 
W.H. Press, Comments Astrophys. 7, 103 (1978); 
Van der Ziel, Adv. Electronics and Electron Phys. 49 (1979) 225. 

[3] R. Botet and Ploszajczak, , in Complex Fluids, p. 303, edited by L. Garrido, 
Springer Verlag (Berlin), 1993. 

[4] R. Botet and M. Ploszajczak, Phys. Rev. Lett. 69, 3696 (1992); Phys. Lett. 
B312, 30 (1993); in Proc. Workshop on Fragmentation Phenomena, Ecole 
de Physique Theorique, Centre de Physique des Houches, 12 - 17 April 1993, 
GANIL Preprint P 93/14 (to be published); Acta Phys. Pol. B25, 353 (1994). 

[5] The third class, associated with the oo - cluster phase, cannot be described in 
terms of the scale-invariance. 

[6] R. Botet and M. Ploszajczak, "New trends in kinetic fragmentation theory", J. 
Mod. Phys. E (in print), GANIL Preprint 94-14, 

[7] An absence of the size dependence in pp is equivalent to the time-independence 
of the fragmentation probability, which in turn follows from the criticality of 
the general FIB model at the shattering transition P, |[] . 

11 



[8] S.P. Obukhov, in Random Fluctuations and Pattern Growth: Experiments and 
Models, edited by H.E. Stanley and N. Ostrowsky, NATO Advanced Study 
Institutes, Ser. C, Vol. 157 (Kluwer, Dordrecht, 1988). 

[9] H.J. Jensen, K. Christensen and H.C. Fogedby, Phys. Rev. B40, 7425 (1989). 

[10] K. Christensen, H.C. Fogedby and H.J. Jensen, J. Stat. Phys. 63, 653 (1991). 

[11] J. Kertesz and L.B. Kiss, J. Phys. A23, L433 (1990). 

[12] Z. Olami, H.J.S. Feder and K. Christensen, Phys. Rev. Lett. 68, 1244 (1992). 
[13] M.J. Jensen, Phys. Rev. Lett. 64 (1990) 3103. 



12 



Figure captions 



Fig. 1 

In the upper part, the normalized average size of the cluster a is plotted as a func- 
tion of time for the FIB automaton with iV = 512 and pp = 0.1 . The solid line 
shows the asymptotic value a = 4/9 when t — > oo . The bottom part of the fig- 
ure shows the cluster-size distribution of the FIB automaton in the steady state for 
Pp = 0.1 (filled circles) and 0.25 (filled triangles) in the SOC phase, for pp = 0.5 
(open circles) at the critical point of the branching process and for pp = 0.75 (open 
squares) above it. The slope of the power-law distributions n(k) is r = 2pp , in 
accordance with results found earlier in the FIB model Q . Notice a total absence 
of the finite-size effects even for k = N . The data correspond to 10 6 fragmentation 
avalanches with N = 256 . 

Fig. 2 

Lifetime distributions of the fragmentation avalanches for different values of pp . 
Significant deviations from the power-law distribution for large times are seen for 
Pp > 0.5 and are due to the appearance of 'infinite' branching trees which reach the 
cutoff for monomers. The data correspond to 10 6 fragmentation avalanches with 
iV = 2048 . The same signatures are used as in the bottom part of Fig. 1. 

Fig. 3 

Power spectra corresponding to the lifetime distributions shown in Fig. 2. The same 
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signatures are used as in the bottom part of Fig. 1. 



Table 1 

Exponents for different values of the fragmentation probability pp ■ In the last col- 
umn, the sum of calculated exponents (3 and r is given. (3 + r — 4 for sandpile 
cellular automata as verified by Christensen et al[|Kj. 



Pf 


P 


a = t 


P + T 


0.05 


1.0 ±0.1 


3.0 ±0.3 


4.0 ±0.4 


0.10 


1.16 ±0.05 


2.81 ±0.05 


3.97 ±0.1 


0.25 


1.50 ±0.03 


2.53 ±0.10 


4.03 ±0.1 


0.40 


1.91 ±0.05 


1.90 ±0.10 


3.81 ±0.1 


0.50 


2.00 ±0.02 






0.75 


2.00 ±0.05 
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